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ABSTRACT 

This paper employs the Navier-Stokes equations, which are modified to include a rotational non- equilibrium 
relaxation model to analyse the heat transfer, drag, and shock stand of distance for hypersonic flow past an 
axisymmetric blunt body for various levels of rarefaction including the rotational non-equilibrium effect. The 
commercial flow solver FLUENT is used to calculate the numerical solutions for laminar viscous hypersonic flow 
past a blunt body at Knudsen number (k t2 ) in continuum-transition regime with and without rotational non¬ 
equilibrium. The effects of rarefaction in the continuum-transition regime are modelled by applying the Maxwell 
velocity slip and temperature jump boundary conditions on the surface. The effects of the rotational non¬ 
equilibrium terms are discussed in this paper for both the continuum K L -, = 0 .land no slip flow regimeK n < 0.1. 
Study is performed for (M m = 6.5, 7, 7.1 and 8). Finally we conclude by comparing drag and heat transfer changes 
by changing K n number. Effect of increase in Mach number is also studied and concluded. 
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1. INTRODUCTION 

During the first stage of atmospheric re-entry, vehicles undergo high heat loads and increased drag due to the high 
velocity and temperature. At these altitudes, the flow can no longer be considered to be in the continuum regime 
due to the effect of reduced atmospheric density, referred to as rarefaction. Many commercial flow solvers are not 
able to accurately predict these flows outside the continuum regime, which challenges the design and optimization 
of space vehicle geometries under those conditions. Simulating these flow conditions is difficult; however, creating 
a flow solver capable of accurately predicting these flows are much more cost effective compared to wind tunnel 
testing in an actual hypersonic environment. 

Rarefaction is characterized by the Knudsen number, a non-dimensional parameter given by: 

x * = i (1) 

Where A F ?J the Knudsen number isA is mean free path of the molecules in free stream and L is the 
characteristic length. 
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The flow is considered to be within continuum regime when Knudsen number is very small <0.002. An accurate 
solution can be calculated in continuum regime using Navier-Stokes equation without slip conditions. These equations 
break down as rarefaction effect increase. It is possible to modify NS equations in order to obtain accurate results in the 
continuum-transition slip regime. This is defined as 0.002<K n <0.1. This is done by implementing Maxwell Slip wall 
boundary conditions. Another significant effect that needs to be considered when hypersonic vehicles re-enter is the 
atmosphere effect of a rotational non-equilibrium. 

In this study, ANSYS Fluent flow solver is used to implement the Maxwell a slip boundary conditions as well as 
the rotational thermal non-equilibrium for the calculation of rare flow fields. The effects of a two-temperature thermal non¬ 
equilibrium model are analysed for blunt leading edge geometry. An axisymmetric blunt body at various Knudsen 
numbers. The geometry is then optimized for the reduction of heat transfer and drag at a chosen Knudsen number with the 
temperature model. 

2. BACKGROUND 


Hypersonic flows have the properties of Shock stand-off distance, Real gas, entropy change, Low density, thermal non¬ 
equilibrium and molecular dissociation. The density of gas behind shock wave increases as Mach number increases. This causes 
the shock wave stand-off distance to decrease with Mach number. This decreased shock stand-off distance mixes with high 
temperature due to kinetic energy causes internal energy. This heat transfer can be reduced by reducing the leading edge radius 
and the drag on the body also increases. We can reduce both drag and heat transfer with the optimization of the body. 

The high temperature associated with hypersonic flows can also increase thermal non equilibrium in rotational 
and vibration modes. For a diatomic gas such as nitrogen, the vibration mode activates at 3300K and become significant 
above 5000K temperature. 

The 3D navier stokes equation for a diatomic gas in conservation law form in Cartesian coordinates. They can be 
written as 

SQ + ^ + Mfc ^=5 

dt &x dy dx dx ds (2) 


Tp a - 


The stress tensor, Ljl and heat flux vector ^ Ljl are obtained from Chapman-Enskog expansion. 

T i; = t J + rf-j + tJ + ■■ ■ + tJ + 0 (JQ+ l ) 

Qij = Qij + Qij + q? } - 4—+ o(A r " +1 ) 


(3) 


(4) 


The navier stokes equations are obtained by taking first two terms of expansion, second order terms, in this stress 
tensor and heat flux vector expansion becomes 


T ..= t p.+ x i = - lot * 

c -> ^ & Xi 


0,1 k9T 

q tJ =^ + q tI = -T- 


(5) 


( 6 ) 


The local viscosity of the flow is calculated using Sutherland’s law 
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H = Pref 



(7) 


H. r _ r T ra f 

Where ' ' ~ J and 'are the reference viscosity and temperature; S is the Sutherland temperature of the gas; 
and T is the local temperature. In order to close the system of equations an equation of state is used, in this case, the 
ideal-gas law. 


3. RAREFACTION EFFECTS 


As space vehicles travel through the atmosphere and into the lower earth orbit, the density of the atmosphere decreases. 
This reduced density is called the rarefaction. The degree of rarefaction is characterized by the Knudsen number: 


K 


n 


A 

L 


( 8 ) 


K A. 

Where 71 the Knudsen number is is mean free path of the molecules in free stream and L is the characteristic 


length. 


The mean free path is the average distance travelled by a particle in the flow between collisions. In general, the 
mean free path of a molecule is given by Equation 
JRT 

(9) 


^ = — 




Where R is the universal gas constant, T is the free stream temperature, d is the average molecule diameter, NA is 
Avogadro's number and P is the free stream pressure. 

In the case of standard air, Equation 9 can be simplified to Equation 10. 



Where ^ is the free stream density and ^ is the viscosity. 

It can be seen in Equation-10 that as the free stream density decreases, the mean free path increases and thus 
increasing the effect of rarefaction. At extremely flow Knudsen numbers (Kn< 0.01), the flow is in the continuum 
regime. Some equations used to model flow, such as the Euler and Navier-Stokes equations, are valid only in the 
continuum regime. As the Knudsen number increases, flow becomes rare field and enters the continuum-transition 
regime (Kn <10) and the Navier-Stokes equations are no longer able to accurately predict the flow. The validity of 
Navier-Stokes equations can be extended to Kn <0.1 by using Maxwell's slip flow boundary conditions [5] on the 
wall. The mathematical equations valid for each flow regime are shown in figure 1.1. 
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Figure 1: Flow Regimes and Appropriate Equations for Modelling the Flow at a Given Knudsen Number. 

4. DESIGN AND MODELLING 

In this chapter, the axisymmetric blunt body simulations in hypersonic flow are performed for four Knudsen numbers 
(0.002, 0.01, 0.05 and 0.1) at (Af M = 6.5, 7.1, 7.5 & 8). The Maxwell's slip boundary conditions for slip velocity and 
temperature jump at the wall are employed in the simulations. The maximum heat transfer and the drag coefficients are 
obtained for each simulation. The geometry and free stream flow conditions used in the simulations correspond to the 
Lobby’s ballistic range experiments and are as given in Ref. [5]. 


The leading edge radius of the sphere is 6.35 mm. Since the body shape is axisymmetric, only a quarter of the 
sphere is used for simulation of the flow. The quarter geometry of the sphere used is shown in figure 2. In the simulations, 
a cylindrical body shape is attached to the rear face of the sphere also shown in figure 2 to exclude the wake effects. A 
structured mesh around the geometry along with a cylindrical body attached is generated using ICEM-CFD. A slice of the 
final mesh is shown in figure 2. A mesh was chosen based on an extensive mesh validation study conducted by Seager [3], 
who employed several types of meshes of varying densities shown in figure 3 and concluded that the segmented mesh (35- 
105-10)x400 was the most efficient mesh with acceptable accuracy. 



Figure 2: Axis Symmetric Blunt Geometry and Mesh. 
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Figure 3: Validation Study Conducted. 


5. NUMERICAL STUDY 

Second order implicit algorithm used for spatial accuracy and dual time stepping transient flow field with multi-grid 
convergence acceleration is employed and the transient flow field is converged at every global time step. The fluid is 
treated as compressible perfect gas with composition of standard air. The k — £ shear-stress transport (SST) model is used 
for turbulence modelling. The body surfaces of the inlet are assumed to be no-slip and adiabatic. The free stream flow 
conditions are obtained from Lobby’s experiment and are given in Ref. [5]. The free stream Mach numbers are 6.5, 7.1, 7.5 
and 8; the temperature is 293K; the air species distribution is 23.3 percent 0 2 and 76.7 percent N 2 . The wall is kept at a 
constant temperature of 1000 K. Since the species distribution is that of standard air, the mean free path given by Equation 
10 is valid for calculating the Knudsen number. The mean free path is calculated from the Knudsen number and the 
characteristic length of the body based on the diameter of 0.0127 m. The free stream density is calculated using Equation 
10 and the free stream pressure is obtained using the ideal gas law. The free stream pressures required to get the desired 
Knudsen numbers are given in Table 1. Residuals are continuously monitored for continuity, k, axial-velocity, and radial- 
velocity. The convergence criterion is that all of these residuals are dropped below 10‘ 3 with the mass flow rate and the 
Mach number of mass-weighted average at throat of inlet retaining constant. 


Table 1: Free Stream Parameters for Blunt Body 


Knudsen Number 

Pressure (Pa) 

Mean Free Path (m) 

Characteristic Length 

0.002 

245.59 

2.54E-5 

0.0127 

0.01 

49.854 

1.27E-4 

0.0127 

0.05 

9.824 

6.35E-4 

0.0127 

0.1 

4.912 

1.27E-3 

0.127 


6. RESULTS AND DISCUSSIONS 

6.1. Pressure Study 

In this study, we observed, as the Mach number increases the pressure also increases, however, an increase in k n number 
resulted in decreased density of atmosphere which decreased pressure as well. We noticed that for Mach = 6.5 at Knudsen 
number 0.1 has the lowest pressure which means too less drag. When it comes to K n number 0.002 pressure increased 
parallel drag also increased. Similarly the rest Mach numbers also vary accordingly by change in K n number. 

Pressure obtained for Mach Number=6:5 

The change in Mach number hasn’t have an effect on Density of the Fluid; so the inputs were given which were calculated 
for increasing in k n value. So using this for k n = 0.002, the pressure generated was around 1117.907KN which is the peak 
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pressure recorded. Later on, it is decreased by increasing k n = 0:01; the drag generated was 2274.598KN. It was further 
diminished to 4480.643KN by increasing the value of k n = 0:05. Finally for k n = 0.1, the pressure generate was 
2239.041KN, which is too lower compared to all the above. 



Figure 4: Pressure Variations for Different Knudsen Number for Mach-6.5. 



Figure 5: Pressure Variations for Different Knudsen Number for Mach-7.1. 

Pressure obtained for Mach number = 7:5 

For Mach number M=7.5 and K n = 0:002, the pressure was recorded around 1588.394KN and it was decreased with the 
increase of the K n = 0:01; the pressure generated was 323189:5N. Further, the pressure was diminished to 63664:29N by 
increasing the K n = 0:05. Finally, the process was stopped by increasing the K n = 0:1 which generated pressure of 
31817:87N and concluded by less pressure observations. 



Figure 6: Pressure Variations for Different Knudsen Number for Mach-7.5. 


Pressure obtained for Mach number = 8:0 

With the increase in pressure generation for each K n value, there is an effect of increase in Mach number; Where for K n 
= 0:002, the pressure generated was about 2411835N and later by increasing the K n value, the pressure generation was 
decreased to 490734.6N for K n = 0:01. Further, pressure was decreased to 96669:16N for K n = 0:05 which is as low as 
possible. Even though the pressure generated was low, it further diminished to K n = 0:1. 
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Figure 7: Pressure Variations for Different Knudsen Number for Mach-8. 


6.2. Temperature Study 

Temperature plots shows decrease in temperature with the increased Knudsen number which results the shock stand-off 
distance will reduce with the increase in Knudsen number results to decrease in drag during flight. Temperature increased 
with the increased Mach number. We observed that temperature decreased till M=7.1 and later on temperature rose with an 
interval of 200K. This shows the heat flux generation on the surface of the vehicle. 

Temperature Variation for M = 6:5 

During hypersonic speed, the heat flux generation was the major drawback which is generated around 2781.378K for K n = 
0.002, which is the lowest heat flux generated, and later on, it increased to 2810.784K for K L -,=0.01. It was observed that 
the temperature was further increased to 2931.116K and finally concluded with 3931.32IK, which the peak temperature is 
obtained on the blunt body. As the density of the fluid increases, the spacing between molecules increases, which leads to 
high collaboration of the molecules. This bonding, when a molecule travels and touches to another molecule with high 
velocity, the temperature increases drastically. Due to this, increase in temperature was recorded for all the K n numbers. 



Figure 8: Temperature Variations for Different Knudsen Number for Mach-6.5. 


Temperature variation for M = 7.1 

As the speed of the vehicle increases, the molecule collaboration also increases and due to loosely packed molecules, the 
velocity distribution of molecules is higher. This increases the temperature to the peak with the increase in K n value. We 
recorded 3260.558K for K n =0.002 and it increased to 3299.171K by increasing K n =0.01. The temperature was further 
increased to 3447.99K for K n =0.05 and finally ended with peak temperature of 3544.279K. 



Figure 9: Temperature variations for different Knudsen Number for Mach-7.1. 
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Temperature variation for M = 7.5 

As stated above, temperature increases with an increase in K n value. For M = 7:5, we observed lowest temperature of 
3608:079K for K a = 0:002 and it elevated slightly by increasing K n to K n = 0:01; the temperature obtained was 
3647:396;. Later on K^value was increased to K n = 0:05 also increased the temperature to 3818:46K and finally ended by 
reaching peak temperature to 3818:46K for K n = 0:1 



Figure 10: Temperature Variations for Different Knudsen Number for Mach-7.5. 


Temperature variation for M = 8:0 

Similarly the temperature rises due to velocity increase by the vehicle. It is because of loosely packed molecules, and these 
molecules collaborate each other with high speed. Our study resulted with a recording of 4311:517N for K a = 0:002 and it 
was increased to 4112:455N for K n = 0:01 and later on it was further increased to 4061:8N by increasing K n =0.05. 
Finally, it was concluded by recording peak pressure of for = 0:1. 



Figure 11: Temperature Variations for Different Knudsen Number for Mach-8. 


7. CONCLUSIONS 

The study on Drag given in below table gives us decrease in drag with an increase in Knudsen number. But, as the Mach 
number increases, drag also increases from the table below, it is inferred that there is a drastic reduction of drag for K n = 
0.1, but when it comes to K n =0.002, we observe only nominal decrease. 

Table 2: Drag Estimation for Various Knudsen Number Increasing with an Increased Mach Number 


Mach number 

K n =0.002 

K n =0.01 

K n =0.05 

© 

II 

a 

X 

6.5 

1.6108603 

0.40040256 

0.11345538 

2.4E-4 

7.1 

1.9064469 

0.46519763 

0.13346742 

2.9E-4 

7.5 

2.1176512 

0.51583894 

0.14780244 

3.3E-4 

8 

2.3983547 

0.58319674 

0.16767936 

3.7E-4 


Pressure plots of various Mach numbers indicate reduction in pressure for increase in Knudsen number and we 
also observed that drag reduced with an increased Knudsen number. However, pressure increases with an increase in Mach 
number, but figure 16 shows that pressure ends up at same point for all the Mach numbers. 
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Figure 12: Pressure Plot for M=6.5. 



Figure 13: Pressure Plot for M=7.1. 
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Figure 14: Pressure Plot for M=7.5. 



Figure 15: Pressure Plot for M=8 and Negligible. 
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Figure 16: Mach Number vs K„=245. 

In the figure-16, an observation shows that peak drag is generated for low Knudsen number. Hence the drag is 
reduced with an increase in Knudsen number. Drag of 2:39835N for initial condition of K n = 245; later on K n value is 
increased from 245to k u = 49:8 which resulted with 0:58319N; further study is performed for K n = 9:8 which gave 
0:16767N. The process was closed after K n = 4:9 which obtained a drag of 0:00033N which is too small. 
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